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Abstract. The the initial stage of a relativistic heavy ion collision can be described 
by a classical color field configuration known as the Glasma. The production of quark 
pairs from this background field is then computed nonperturbatively by numerically 
solving the Dirac equation in the classical background. The result seems to point 
towards an early chemical equilibration of the plasma. 



PACS numbers: 24.85. +p, 25.75.-q, 12.38.Mh 
1. Introduction 

At large energies (small x) the hadron/nucleus wavefunction is characterized by a 
saturation scale Q s arising from the strong nonlinear interactions between the wee 
partons, mostly gluons. If the energy is high enough so that Q s ^> Aqcd, weak coupling 
methods can be used to decribe physics at transverse momenta ~ Q s . In the context of 
ultrarelativistic heavy ion collisions this means that one can hope to understand not only 
hard probes, but the bulk of particle production in terms of weak coupling, deconfmed 
physics. Note that although the coupling is weak, the saturated color fields are strong 
~ 1/g, and the physics must still be treated nonperturbatively. 
The purpose of this talk is to explore the initial state of a heavy ion collision in this 
framework. We will first describe what has been called the Glasma PP, a highly coherent 
classical field in the initial stage of the collision. We will then move on to study how 
quark pairs can be produced from this, to a first approximation, purely gluonic system, 
moving it closer to a chemically equilibrated plasma of both quarks and gluons [21 0] • 



2. Glass and Glasma 



The small x wavefunction of a hadron or nucleus, characterized by nonperturbatively 
large color fields, can be described in terms of a classical Weizsacker- Williams (WW) 
field radiated by the hard, large x, sources. Because of their high speed and Lorentz 
time dilation, the hard degrees of freedom are seen by the low x fields as slowly evolving 
in lightcone time. They can therefore be thought of as classical, static (in light cone 
time) sources for the small x fields jjj. This effective description has been called the 
Color Glass Condensate. The earliest stage of an ultrarelativistic heavy ion collision 
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is a coherent, classical field configuration of two colliding sheets of Color Glass. This 
first fraction of a fermi of the collision, in transition from two sheets of Color Glass into 
eventually a quark gluon plasma, is what we refer to as the Glasma pQ. 

The color fields of the two nuclei are transverse electric and magnetic fields on the 
light cone. The glasma fields left over in the region between the two nuclei after the 
collision at times 1 < r < 1/Q S are, however, longitudinal along the beam axis. One 
way of understanding these field configurations is the following. Let us work in light 
cone gauge, so that each nucleus, when going past a point on the beam axis with no 
gauge field before the collision, leaves behind it a pure gauge field (see Fig. [TJ. One can 
define an effective chromoelectric and chromomagnetic charge density by separating the 
nonlinear parts of the vacuum Gauss law and Bianchi identities 

[A, E?] = and [A, B l ] = (1) 

as 

diE* = Pc = ig[A\ E 1 } and = p m = ig[A\ B% (2) 

Now we can interpret the interaction of the WW chromoelectric and -magnetic fields 
of the nucleus on the x + -light cone with the pure gauge field left behind by the other 
nucleus as an effective chromoelectric and -magnetic charge density left behind on the 
light cone. An exactly opposite charge density is left behind on the other sheet, leading 
to a longitudinal chromoelectric and -magnetic field between the sheets \. This structure 
is illustrated in Fig. |2j 

As the initial condition has a chromoelectric and -magnetic field both in 
the longitudinal direction, it also has a nonzero Chern-Simons charge density ~ 
g 2 e M " ,pa tr F^Fpcr, fluctuating on transverse length scales of the inverse saturation 

f Note that the initial fields being longitudinal in along the beam axis direction is in no contradiction 
with the lowest order perturbative description of the process as gg — ► g scattering, because the 
longitudinal (with respect to the beam axis) fields are perpendicular to the momentum of the gluon 
being produced. The initial polarization state of this gluon is, however, a very particular one. 
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Figure 2. The WW fields of the two nuclei before and after the collision. Before 
the collision there are only transverse fields on the sheets. After the collision the 
interaction of these fields with the pure gauge field of the other nucleus leaves behind 
an effective electric and magnetic charge density (the dots on the figure) on the sheet, 
and a longitudinal electric and magnetic field between these effective charges. 



momentum. This can naturally lead to large P- and CP-violations in these domains, 
which could have interesting consequences as parity-odd global observables in the 
produced mesons El HI IE] • If the gauge field configurations are exactly boost invariant, 
no topologically nontrivial configurations (sphaleron transitions) are allowed [0], and 
the effect is suppressed. But as giving up the assumption of boost invariance in the 
numerical calculation is becoming more feasible PHI one will be able to understand 
this phenomenon better. 

How, then, does the glasma decay into plasma, i.e. how can we move from a 
description of the initial state in terms of classical fields into one formulated for quantum 
particles? The glasma fields depend on the transverse coordinate on a length scale 
of order 1/Q S . Thus to lowest order (in the coupling a s or in h) the fields simply 
radiate away as gluons with p T ~ Q s . As the system expands the fields are diluted 
and can be treated as particles. This lowest order production is the contribution that 
is computed in the numerical computations of gluon production in heavy ion collisions 

da nam en EE] • 

In the rest of this paper we will be concerned with the next order in g or H, 
namely quantum pair production from the classical background field. Both gluons and 
quarks carry color charge and can be produced in pairs from the background field. The 
production of gluon pairs is the first quantum correction to the classical background field. 
To consistently compute it one must take care to avoid double counting the gluons that 
are already effectively included in the classical background field. A calculation of gluon 
pair production must be done consistently with the high energy renormalization group 
evolution of the sources. This computation has not yet been precisely formulated. 
On the other hand, computing quark pair production from the classical background 
field is more directly relevant for the phenomenology of heavy ion collisions and easier 
to formulate, if not necessarily trivial to carry out nonperturbatively, because it only 
involves the lowest order classical gluon field. 
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3. Background field from the MV model 

The hard degrees of freedom are modeled as classical sources on the light cones: 

r = ^ + P(1) (x T )5(x-) + 5^ P{2) (x T )S(x + ). (3) 

The Weizsacker- Williams fields describing the softer degrees of freedom can then be 
computed from the classical Yang-Mills equation [D^, F^ u ] — J u . In the light cone 
gauge the field of one nucleus is a pure gauge outside the light cone (see Fig. HJ) 

4 2) = -U^dM 2) , with t/ (m) (x T ) = exp (4) 

In the original MV model, which we shall be using here, the color charge densities are 
stochastic Gaussian random variables on the transverse plane 

<p a (x T )/(y T )> = <? 2 /iW(x r - y T ), (5) 

where the density of color charges g 2 fi is, up to a numerical constant and a logarithmic 
uncertainty, proportional to the saturation scale Q s . 

The initial conditions for the fields in the future light cone between the two colliding 
sheets were derived and the equations of motion solved to lowest order in the fields in 
Refs. [T7J HH] (see also Ref. [TH] for the same calculation in covariant gauge and Ref. [20] 
for another formulation of the same lowest order result.) This initial condition has a 
simple expression in terms of the pure gauge fields (j3J) of the two colliding nuclei: 

^ = 40 +4) and ^ = f[4)>4d- ( 6 ) 

The solution of the Yang Mills equations beyond leading order is not known 
analytically. They can, however, be solved numerically on the lattice. The numerical 
setup was formulated in Ref. [12] and used to calculate the energy density and gluon 
multiplicity corresponding to the glasma fields in e.g. Refs. [121 EB1 UH El EE]- It is 
these numerically computed color field configurations that we will use as the background 
field for computing quark production in the following. 

4. Pair production 

Let us then move on to the computation of quark pair production from the background 
color field described above starting with a few trivial remarks: 

• We are assuming such a high collision energy that the nuclear wavefunction is 
completely gluonic and the process gg — > qq dominates over quark production from 
the sea and valence distributions by processes like qg — > qg. This might not yet 
be a very good approximation at RHIC energies [21] • We are also neglecting the 
backreaction of the produced pairs on the gluon fields. 

• When the computation of gluon production to the lowest order is done 
nonperturbatively to all orders in the strong field the multiplicity and energy density 
of the gluons is infrared finite. But for the computation of quark pair production 
it is not a priori evident whether the limit m q — > is finite. 



Chemical composition of the decaying glasma 



5 




x 



X 





1 I 1 1 1 1 1 




: 






q = (0,2) 






q = (4,4) 






-- - q=(10,5) 






q = (20,13) 










s 

* 


N 

s 







Ay 



Figure 3. Left: Spacetime in Dirac equation. In the case when the nuclei move 
exactly at the speed of light the solution decomposes into two distinct paths depending 
on whether the x + or x~ light cone is encoutered first. In the Abelian case these 
correspond to the t and u channel terms of the amplitude. Right: amplitude as a 
function of rapidity y p — y q for different transverse momentum modes qy. 



• If quark production is dominated by saturation physics, one could expect production 
to be flavor blind for m q <C Q s . Thus strangeness would, from the beginning, be 
equilibrated with the light flavors. 

Our method of computing quark pair production by solving the Dirac equation is 
explained in more detail in Ref. [2]. One starts, at t — ■> — oo with a negative energy plane 
wave solution of the free Dirac equation, e iq x v(q) (see Fig. EJ). The Dirac equation is 
solved forward in time until t — > oo, and the solution projected on positive energy plane 
wave states e~ ip ' x u{p). The amplitude is then interpreted as the amplitude to produce 
a quark pair from the classical background field. Solving the Dirac equation forward in 
time is equivalent to computing the full retarded (not Feynman) Dirac propagator in 
the background field. This means [22] that one is computing the expectation value of 
the number of pairs produced, not the cross section to produce exactly one pair (which 
would be related to the Feynman propagator). The formal derivation of this procedure 
is given in Ref. [22], but the physical interpretation is perhaps more intuitive in terms 
of the "Dirac sea". Picturing the vacuum as a Dirac sea with all the negative energy 
states filled, one is taking a component of the "in" vacuum (at t — > — oo) and computing 
its overlap with an "out" state (at t — > oo) of a quark with momentum p. If this overlap 
is nonzero, the background field has lifted a negative energy solution from the Dirac 
sea to positive energy, leaving a hole in the sea. This process is then interpreted as the 
creation of a pair of fermions. 

Because the background field of one nucleus is a pure gauge (see Figs. [T] and 0J), the 
Dirac equation can be solved analytically up to the future light cone (x + = 0, x~ > 
and x~ = 0, x + > 0). In the Abelian case also the field inside the future light cone 
is a pure gauge, and the whole computation (electron positron pairs in ultraperipheral 
collisions of two heavy ions) can be carried through analytically j2S] . In our case the field 
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inside the future light cone is only known numerically, and so we must solve the Dirac 
equation numerically starting from an analytically derived initial condition at r = 0. It 
turns out that this is most conveniently done in the (r, z, XT)-coordinate system. The 
natural choice for the temporal coordinate is r, because the initial condition for the 
numerical computation is given at r = 0, and thus no other time coordinate would lead 
to a pure initial value problem. For the longitudinal coordinate one must (unlike the 
three dimensional computation of the color field in Refs. fUl^J) choose a dimensional 
coordinate in order to represent the longitudinal momenta (q^ on a coordinate space 
lattice at r = 0. 

It was shown in Ref. that in the lowest this calculation of pair production 
reduces to a standard expression in kr-factorized perturbation theory. The "pA" 
case, when only one of the classical color fields is treated to lowest order, can also 
be solved analytically |2*H] . A kr-factorized perturbative formalism has also been used 
to compute production of heavy quarks in from the "Color Glass Condensate" gluon 
distributions e.g. in Ref. |2E]- A way of approaching the problem from the other 
limit is presented in Ref. j2Z], where, using the WKB approximation, pair production 
is computed nonperturbatively in a short longitudinal color field pulse, neglecting 
the magnetic field. Because this computation also neglects the transverse coordinate 
dependence of the field, it cannot reduce to the same perturbative expression in the 
weak field limit. 

Let us then finally move to the results of the numerical computation [3]. The 
computation has three independent numerical parameters, the color charge density g 2 n 
that determines the strength of the background field, the nuclear radius Ra giving the 
size of the system in the transverse plane, and m, the quark mass. Figure El shows 
the amplitude, integrated over the momentum of the quark but for different transverse 
momenta of the antiquark, as a function of the rapidity difference Ay = y p — y q between 
the quark and the antiquark. Because the background field is boost invariant, this 
amplitude only depends on Ay and not the rapidity of the whole system (y p + y q )/2. 
When also integrating over Ay one gets the spectrum as a function of the transverse 
momentum of the antiquark, shown in Fig. 0] for different strengths of the background 
color fields. On a finite lattice the large part of the spectrum is affected by the 
proximity of the lattice cutoff ~ 1/a, which is demonstrated in Fig. El where the same 
spectrum is plotted for different values of the lattice cutoff. Although physically the 
quarks can only be interpreted as on shell particles after a formation time ~ 1/mx, the 
projection on positive energy estates can be done at arbitrarily early times. A shown 
in Fig. El the amplitude is close to its final value for very early times close to the light 
cone. 

According to conventional wisdom the initial state of a heavy ion collision is 
dominated by gluons. Assuming that the subsequent evolution of the system conserves 
entropy this would mean ~ 1000 gluons in a unit of rapidity. In the classical field 
model this corresponds [IB| to g 2 [i = 2 GeV. Our result seems to point to a rather large 
number of quark pairs present already in the initial state. One could envisage a scenario 




Figure 4. Left: Antiquark transverse momentum spectrum for different values of the 
classical color charge density g 2 /i for m = 0.3 GeV. The straight lines to (mis)guidc the 
eye are ~ e ~ q - T ^ 9 W. Right: The antiquark spectrum as a function of the transverse 
mass \J qr 2 + m 2 for different masses and g 2 jx = 1 GeV. 
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Figure 5. Left: Amplitude with m = 0.3 GeV and g 2 /j, = 1 GeV for different 
lattice sizes, corresponding to different values of the ultraviolet lattice cutoff. Right: 
Amplitude as a function of the time at which the projection to positive energy states 
is done. The first four data sets are for g 2 /j, = 2 GeV and the last one (marked with 
an asterisk) for g 2 /i = 1 GeV. 



where, for g 2 fi = 1.3 GeV, these 1000 particles could consist of 400 gluons, 300 quarks 
and 300 antiquaries (take the lowest curve from Fig. Eland mutiply by N{ = 3). This 
would be close to the thermal ratio of N q /N g = 9N { /32. 

5. Conclusions 

We have pointed out some known, but not always fully appreciated, features of the 
classical field description of the early stages of a heavy ion collision. The "Glasma" fields 
are initially longitudinal chromoelectric and chromomagnetic fields, of equal magnitude. 
These fields can be thought of as arising from effective chromoelecric and -magnetic 
charge densities on the transverse planes caused by the color rotation of the Weizsacker- 
Williams fields of one nucleus in the pure gauge field of the other. 

We have then studied quark pair production from classical background field of 
McLerran-Venugopalan model studied by solving the 3+1-dimensional Dirac equation 
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numerically in this classical background field. We find that the number of quarks 
produced large, potentially leading to a very early chemical equilibration of the quark 
and gluon degrees of freedom. But in order to put this conclusion on a more theoretically 
sound footing one must also compute gluon production to the same order in the 
couopling a s , which has not yet been done. Our numerical method does not yet allow 
us to extend our method to large quark masses or transverse momenta. 
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